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X-RAY MICROLENSING IN RXJ1131-1231 AND HE1104-1805 
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ABSTRACT 

We present results from a monitoring campaign performed with the Chandra X-ray Observatory of 
the gravitationally lensed quasars RX J1131 — 1231 and HE 1104—1805. We detect significant X-ray 
variability in all images of both quasars. The flux variability detected in image A of RX J1131— 1231 
is of particular interest because of its high amplitude (a factor of ~ 20). We interpret it as arising 
from microlensing since the variability is uncorrelated with that of the other images and the X-ray flux 
ratios show larger changes than the optical as we would expect for microlensing of the more compact 
X-ray emission regions. The differences between the X-ray and optical flux ratios of HE 1104—1805 
are less dramatic, but there is no significant soft X-ray or dust absorption, implying the presence of 
X-ray microlensing in this system as well. Combining the X-ray data with the optical light curves we 
find that the X-ray emitting region of HE 1104—1805 is compact with a half-light radius <6r s , where 
the gravitational radius is r g = 3.6 x 10 14 cm, thus placing significant constraints on AGN corona 
models. We also find that the microlensing in HE 1104—1805 favors mass models for the lens galaxy 
that are dominated by dark matter. Finally, we better characterize the massive foreground cluster near 
RX J1131 — 1231, set limits on other sources of extended X-ray emission, and limit the fluxes of any 
central odd images to be 30-50 (3cr) times fainter than the observed images. 

Subject headings: galaxies: active — quasars: individual (RX J1131-1231) — quasars: 
individual (HE 1104-1805) — X-rays: galaxies — gravitational lensing 



1. INTRODUCTION 

Resolving the emission regions of distant quasars is be- 
yond the current capabilities of present-day telescopes, as 
direct imaging of accretion disks requires angular resolu- 
tions on the order of tens of nanoarcseconds at z ~ 1. Until 
the spatial resolution of telescopes reaches this limit, we 
will have to rely on indirect methods of mapping the emis- 
sion regions of quasars. Such methods include light-travel 
time arguments, reverberation mapping of the broad line 
region (Blandford & McKee 1982; Peterson 1993, Netzer & 
Peterson 1997), reverberation mapping of the Fe Ka emis- 
sion region (Young & Reynolds 2000), occupation mea- 
surements of the central X-ray source by orbiting Comp- 
ton thick clouds (Risaliti et al, 2007), and microlensing 
of the continuum and line emission regions (e.g., Grieger 
et al. 1988 and 1991; Schneider, Ehlers & Falco 1992 
Gould & Gaudi 1997; Agol & Krolik 1999; Yonehara et al 
1999; Mineshige & Yonehara 1999; Chartas et al. 2002a 
Popovic et al. 2003; Blackburne et al. 2006; Pooley et al 
2006, 2007; Kochanek et al. 2007; Jovanovic et al. 2008; 
Morgan et al. 2008a, 2008b). 

We observe microlensing as uncorrelated variability in 
the brightnesses of the images of a lensed quasar, where 
the amplitude of the variability is determined by the size 
of the emission regions (e.g., Lewis et al. 1998; Popovic & 
Chartas 2005). The largest components, the radio, dust 
and optical line emission regions, should show little or no 
microlensing variability. The thermal continuum emission 
from the disk should show greater variability at shorter 
wavelengths corresponding to smaller disk radii and higher 
temperatures, and this has been observed by Poindexter et 
al. (2008) and Anguita et al. (2008). The non-thermal X- 
ray emission is thought to be dominated by inverse Comp- 



ton scattering of UV photons from the disk (e.g., Reynolds 
& Nowak 2003), some of which is reprocessed into Fe Ka 
line emission (e.g., George & Fabian 1991; Mushotzky et 
al. 1993) but the geometry and scale of these emission 
regions is not well understood. 

The effects of microlensing on the X-ray fluxes of lensed 
images have been reported for many lenses at this point, 
including Q0957+561 (e.g., Chartas et al. 1995), RX 
J0911+0511 (e.g., Morgan et al. 2001; Chartas et al. 
2001), H1413+117 (e.g., Chartas et al. 2004, 2007), 
PG 1115+080 (Pooley et al. 2006; Morgan et al. 2008b), 
and RX J1131-1231 (e.g., Blackburne et al. 2006; 
Kochanek et al. 2007; Dai et al. 2008). The differences 
between the X-ray flux ratios and optical flux ratios can 
be used to estimate the source size after accounting for ab- 
sorption by the interstellar medium of the lens (e.g., Dai & 
Kochanek 2008) and the effects of substructures (satellites) 
in the lens halo (Mao & Schneider 1998). Pooley et al. 
(2007) used observed X-ray and optical flux ratios for ten 
four-image lensed quasars to show that the X-ray source is 
much more strongly microlensed than the optical source, 
indicating that the X-ray emission region is more compact 
than the optical. They also found that their estimates 
of the optical source sizes were significantly larger than 
would be expected based on the observed, magnification- 
corrected optical fluxes. Morgan et al. (2008a) measured 
the sizes of the optical emission regions of 11 quasars using 
microlensing to find that the sizes scaled with black hole 
mass as expected from thin disk theory (Shakura & Sun- 
yaev 1973) and had approximately the sizes expected from 
thin disk theory, but confirmed the Pooley et al. (2007) 
result that these sizes, both microlensing and theoretical, 
are larger than expected given the observed optical flux 
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and the same disk model. 

Accurate microlensing estimates of the sizes of the X- 
ray emission regions requires X-ray light curves rather 
than isolated epochs. First, the absolute magnifications of 
lensed images are not well understood because substruc- 
ture in the lens galaxy can modify the fluxes from the pre- 
dictions of lens models (e.g., Mao & Schneider 1998; Dalai 
& Kochanek 2002). While comparing optical and X-ray 
flux ratios provides estimates for the difference in size be- 
tween the emission regions, the lack of a secure estimate of 
the intrinsic flux ratios makes it difficult to set an absolute 
scale. Second, quasars have intrinsic variability that ap- 
pears in the images with relative time delays, which means 
that instantaneous flux ratios are contaminated by the ef- 
fects of time variability modulated by the delays. Third, 
while the source size ultimately determines the amplitude 
of microlensing variability, the instantaneous values are 
a combination of the source size and the location of the 
source in the complex microlensing magnification patterns. 
Monitoring the variability minimizes these problems and 
should lead to far more accurate estimates of the sizes. 

We have been monitoring roughly 25 lenses in the op- 
tical to measure time delays (e.g., Kochanek et al. 2006) 
and to study quasar structure (e.g., Morgan et al. 2008a; 
Poindexter et al. 2008). From this sample we selected 
a small subset with reasonable X-ray fluxes for monitor- 
ing with the Chandra X-ray Observatory. Given a well- 
sampled optical light curve that can be used to determine 
the size of the optical/UV emission region and act as a mi- 
crolensing variability "reference," we can afford to sample 
the X-ray microlensing variability relatively sparsely. In 
this paper we present the X-ray observations of the lensed 
quasars RX J1131-1231 (z s = 0.658) and HE 1104-1805 
(z s = 2.32). The optical data for RX J1131 1231 and 
HE 1104-1805 are presented in Morgan et al. (2006) and 
Poindexter et al. (2007) respectively, and the optical mi- 
crolensing results are presented in Morgan et al. (2008a) 
and Poindexter et al. (2008). Here we carry out a mi- 
crolensing analysis of HE 1104—1805, while a companion 
paper, Dai et al. (2009), does so for RX J1131-1231. We 
discuss the observations and data analysis in §2, study 
the microlensing of HE 1104-1805 in §3 and summarize 
our results in §4. Throughout this paper we adopt a A- 
dominated cosmology with H = 70 km s _1 Mpc -1 , f^A 
= 0.7, and n M = 0.3. 

2. X-RAY OBSERVATIONS AND DATA ANALYSIS 

We monitored RX J1131-1231 and HE 1104-1805 with 
the Advanced CCD imaging Spectrometer (ACIS; Garmirc 
et al. 2003) on board the Chandra X-ray Observatory 
(hereafter Chandra) using short (~ 5 ksec) snapshot ob- 
servations. We obtained 5 epochs for each lens in 2006 
placing the sources on the back-illuminated S3 chip of 
ACIS. We combined these with archival observations of 
RX J1131-1231 in April 2004 and HE 1104-1805 in 
February 2000. A log of the observations that includes 
observation dates, observed count rates, total exposure 
times, and observation identification numbers is presented 
in Table 1. We analyzed the data using the standard soft- 
ware CIAO 4.0 provided by the CXC. We used standard 
CXC threads to screen the data for status, grade, and time 
intervals of acceptable aspect solution and background lev- 



els. We removed the ± 0"25 spatial randomization applied 
to the event positions by the standard pipeline and in- 
stead used the sub-pixel resolution techniques developed 
by Tsunemi et al. (2001) and Mori et al. (2001) in order 
to improve the image resolution. 

The Chandra spectra of RX J1131 — 1231 and 
HE 1104—1805 were fit with a variety of models employing 
XSPEC version 12 (Arnaud 1996). For all spectral models 
of RX J1131 1231 and HE 1104-1805 we included Galac- 
tic absorption due to neutral gas (Dicky & Lockman 1990) 
with column densities of N^— 3.6 x 10 20 cm~ 2 and 4.6 
x 10 20 cm~ 2 , respectively. We use 90% confidence level 
uncertainties unless otherwise stated. 

2.1. Spatial and Spectral Analysis 

For estimating the X-ray counts of images B, C, and D 
of RX Jl 131 — 1231 we extracted events from circular re- 
gions with radii of 1.5 arcsec slightly off-center from the 
images to reduce contamination from nearby images (see 
panel (c) of Fig. 1). To estimate the X-ray counts of image 
A of RX J1131 — 1231, which lies between images B and C, 
we extracted events within a 0.75 arcsec circular region 
centered on A. Aperture corrections were applied to all 
images to account for counts not included in the extrac- 
tion regions. For estimating the X-ray counts of images A 
and B of HE 1104—1805 we extracted events from circu- 
lar regions centered on the images with radii of 1.5 arcsec. 
The backgrounds for RX J1131-1231 and HE 1104-1805 
were determined by extracting events within an annulus 
centered on the mean location of the images with inner 
and outer radii of 7.5 arcsec and 50 arcsec, respectively. 
Contamination between images was minimized and made 
insignificant by the application of the sub-pixel resolu- 
tion technique and the use of appropriately sized and off- 
center extraction apertures. We also corrected for "pile- 
up," an instrumental effect that occurs when two or more 
X-ray photons strike individual or neighboring CCD pixels 
within one frame time. Neglecting these corrections can 
cause spectral distortion, grade migration and distortion 
of the image PSF. We used the forward spectral-fitting 
tool LYNX (Chartas et al. 2000) to estimate the fraction 
of events lost due to the pile- up effect. 

Table 1 summarizes the observed 0.2-10 kcV band 
source counts. The pile-up corrections needed to correctly 
estimate these counts can be significant. For example the 
counts for images B and C in the 2004 April 12 observation 
of RX Jl 131 — 1231 have pile-up corrections of 37% and 
16% respectively. In their analysis of the 2004 observation 
of RX J1131 1231, Blackburne et al. (2006) did not in- 
clude pile- up corrections, so our present results should be 
more reliable. In particular, pile-up usually leads to energy 
spectra that appear harder, and this likely explains the dif- 
ferences in the X-ray spectral hardness ratios reported by 
Blackburne et al. (2006) and our present results. 

In Figure 1 we show the Lucy-Richardson deconvolved 
images in the 0.2-8 keV bandpass of the Chandra obser- 
vations of RX J1131 1231 and HE 1104-1805. To esti- 
mate the relative X-ray image positions of RX J1131 — 1231 
and HE 1104—1805 we modeled the Chandra images us- 
ing point-spread functions (PSFs) generated by the sim- 
ulation tool MARX (Wise et al 1997). The X-ray event lo- 
cations were binned with a bin-size of 0"0246 to sample 
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the PSF sufficiently (an ACIS pixel subtends 0"491). The 
simulated PSFs were fitted to the Chandra data by min- 
imizing the C-statistic formed between the observed and 
simulated images. In Table 2 we compare the X-ray image 
separations to the observed NICMOS values reported by 
Morgan et al. (2006) for RX J1131 1231 and by Lehar 
et al. (2000) for HE 1104-1805. We conclude that the 
X-ray and optical image positions are consistent given the 
estimated uncertainties. 

The deconvolved images of RX J1131 — 1231 and 
HE 1104—1805 do not show any additional lensed im- 
ages. To obtain quantitative limits on possible additional 
images located at the centers of RX J1131 — 1231 and 
HE 1104-1805 we extracted the 0.2-10 keV counts within 
0.4 arcsec circles centered on the images at the mean lens 
position of each quasar. The backgrounds in the central 
source extraction regions are dominated by the contamina- 
tion from the images - the instrumental and cosmic back- 
grounds in the 0.4 arcsec extraction circles are significantly 
lower than the contamination. To estimate the fractional 
contamination per image we used the 0.2-10 keV counts in 
0.4 arcsec circular apertures placed 1 and 2 arcsec North of 
image B in RX J1131 — 1231 and 1 and 2 arcsec Northwest 
of image A in HE 1104-1805. The detected 0.2-10 keV 
counts in the central 0.4 arcsec apertures are consistent 
with the estimated background and contamination from 
the bright images. 

By combining all the data for each lens we set 3a up- 
per limits on the 0.2-10 keV flux of any central, odd 
image of 1 x 10~ 14 and 4 x 10~ 15 ergs s~ 4 cm~ 2 for 
RX J1131-1231 and HE 1104-1805, respectively. These 
limits corresponded to a 3tr limit on the flux ratio of a 
factor of 50 relative to image C in RX J1131 — 1231, and a 
factor of 30 relative to image A in HE 1104-1805. While 
relatively tight, these limits are not strong enough to con- 
strain the central surface density of the lens (Keeton 2003). 
For example, in the models we use for HE 1104—1805 be- 
low, the expected flux ratio between image A and any 
central image is > 10 3 . 

We performed fits to the individual spectra of the im- 
ages of RX J1131-1231 and HE 1104-1805 using events 
in the 0.4-8 keV energy range with a model that consisted 
of a simple power law modified by Galactic absorption. 
Due to the moderate S/N of the spectra, we performed 
these fits using the Cash statistic which does not require 
binning of the data, although if we instead use x 2 statistics 
we find similar results. The best-fit parameters of these 
fits and the unabsorbed 0.2-2 keV and 2-10 keV fluxes 
are presented in Tables 2 and 3. We found no significant 
variability (within the 90% errors) of the photon indices 
(r) of the spectra of the images of RX Jl 131 — 1231 and 
HE 1104-1805 with the exception of the 2004 observation 
of RX J1131— 1231 that showed significant differences of T 
between images A/B/C with T ~ 1.44 and image D with 
r ~ 1.95. We also considered models that included neu- 
tral absorption at the redshifts of the sources or lenses of 
RX J1131 1231 and HE 1104-1805. We do not detect 
any additional neutral absorption at these redshifts. This 
is consistent with the marginal detection of absorption in 
HE 1104-1805 in our earlier analysis (Dai et al. 2006). 

In Figure 2 we show the Chandra image of the 2004 ob- 
servation of the lensed system RX J1131 — 1231 and the 



surrounding field. To reduce background contamination 
and to enhance possible soft extended X-ray emission we 
filtered the image to include only photons with energies 
ranging between 0.4 and 3.0 keV. The image was binned 
with a bin size of 0.5 arcsec and adaptively smoothed with 
the tool CSMOOTH developed by Ebeling et al. (2000). 
CSMOOTH smooths a two-dimensional image with a cir- 
cular Gaussian kernal of varying radius. Wc confirm the 
extended soft X-ray emission reported by Morgan et al. 
(2006) centered at 11 32 1.6, -12 31 6.5 (J2000) and about 
158 arcsec Northeast of image D of RX J1131-1231. This 
emission comes from a foreground cluster at z = 0.1. Mor- 
gan et al. (2006) also reported a possible detection at a 
low significance level (3 a) of extended emission centered 
33 arcsec Southwest of RX J1131-1231 at 11 31 50.1, 
— 12 32 23 (J2000). We combined all six observations of 
RX Jl 131 — 1231 listed in Table 1, for a total exposure time 
of 33.8 ksec, to better limit this possible source. All the ob- 
servations from 2006 were performed in sub-array mode to 
mitigate pile-up and with the reduced field of view they do 
not include the foreground cluster to the Northeast of the 
lens. Figure 3 shows the adaptively smoothed image in- 
cluding only photons with energies from 0.4 to 3.0 keV. We 
find no significant sources of extended emission near the 
lens beyond the foreground cluster. The 0.2-2. keV 3a up- 
per flux and luminosity limits are 5 x 10 -15 ergs cm~ 2 s~ 4 
and 1.3 x 10 42 ergs s _1 , respectively, assuming a cluster 
with a temperature of 1.5 keV at the lens redshift with an 
abundance of 0.3 solar and an extraction radius of 30 arc- 
sec. 

We fit the emission from the foreground cluster using a 
P model for the cluster brightness profile combined with 
a uniform background of 0.005 events per pixel. Prior to 
performing the fit we binned the image in 1" pixels and 
smoothed this with a Gaussian (a — 3"). The fits were 
performed with the CXC software package SHERPA. We 
find that the cluster center is Aa = 150"2 ± 0"5 East and 
AS = 52"3 ± 0"5 North of image D. The smoothed inten- 
sity distribution is nearly round with an ellipticity of e = 
0.10 ± 0.3. The best-fit values for [3 and the core radius 
of the cluster arc (3 = 0.4 ± 0.2 and r = 4"2±°; 4 (8 kpc), 
respectively. 

We extracted the spectrum of this z = 0.1 cluster from 
a 50 arcsec radius circle centered on the X-ray cluster cen- 
ter. This spectrum was fit with a simple model consisting 
of an emission spectrum from hot diffuse gas based on the 
XSPEC model mckal modified by Galactic absorption. We 
obtain best-fit values for the temperature and metal abun- 
dances of 1.2lg' 2 keV and A = 0.5^0 2' respectively (both 
errors are at the 90% confidence level). The 2-10 keV lu- 
minosity of this cluster of galaxies is 1.7 x 10 42 ergs sec -1 . 
These arc consistent with the Morgan et al. (2006) esti- 
mates and the values for the temperature, (3 and the core 
radius are consistent with the observed correlations of clus- 
ters (e.g. Jones & Forman 1999). 

We also stacked the images of the six observations of 
HE 1104—1805 listed in Table 1 for a total exposure time 
of 71.8 ksec. Figure 4 shows the adaptively smoothed im- 
age of HE 1104—1805 including only photons with ener- 
gies from 0.4 to 3.0 keV. We find no significant sources 
of extended emission. The 0.2-2 keV 3cr upper flux 
and luminosity limits are 3.4 x 10~ 15 ergs cm~ 2 s _1 and 
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5.7 x 10 42 ergs s _1 , respectively, assuming a cluster with a 
temperature of 1.5 keV at the lens redshift with an abun- 
dance of 0.3 solar and an extraction radius of 30 arcsec. 

2.2. Timing Analysis 

The 0.2-10 keV time-delay corrected light curves of the 
images of RX J1131-1231 and HE 1104-1805 are shown 
in Figures 5a and 5b respectively. Significant X-ray flux 
variability is detected in all images. RX Jl 131 — 1231 shows 
a mixture of intrinsic and microlcnsing variability. If we 
normalize the first epoch for image D to match the optical 
light curves, we find that images B, C and D roughly track 
the optical light curves and show similarities in their time 
variability. In Figure 6a we show the flux ratios Fa/Fq, 
F B /F C , and F D /F C in the 0.2-10 keV and i?-bands of 
RX J1131-1231. The flux ratio F B /F C is almost constant 
in both the X-ray and optical bands indicating that the 
observed variability in images B and C is correlated and 
therefore mostly intrinsic in origin. Image A shows a com- 
pletely different behavior. It is a factor of <~ 4 fainter in 
the X-rays than the optical relative to C/D in 2004, com- 
parable to the optical in the spring of 2006, and a factor 
of ~ 5 brighter by the fall of 2006. Since the 2006 observa- 
tions come in clusters with temporal separations of order 
the A/B/C time delays, we can be confident that these 
differences are not due to intrinsic variability modulated 
by the time delay. The significant rise of the Fa/Fq ratio 
in the X-ray band compared to the one measured in the R- 
band supports our microlcnsing interpretation of the flux 
enhancement in image A of RX J1131 — 1231. We note that 
the flux ratios reported in our analysis for the 2004 obser- 
vation of RX J1131 — 1231 are different from those reported 
by Blackburne et al. (2006) because of the correction for 
the loss of counts due to pile-up and the correction for the 
spectral distortion due to pile-up. We see no strong evi- 
dence for variations in the spectral hardness, but this is a 
weak statement given the limited count rates. 

HE 1104— 1805is also affected by intrinsic X-ray vari- 
ability. In particular, the peak in the optical flux near 
4100 days seems likely to be associated with a peak in 
the X-ray flux. In Figure 6b we show the flux ratios 
Fa/Fb in the 0.2-10 keV and i?-bands of HE 1104-1805. 
The February/March and October/November epochs are 
roughly separated by the time delay, so we can reasonably 
conclude that the X-ray flux ratio near 4050 days (ob- 
served A) is approximately Fa/Fb— 1 and very different 
from the optical flux ratio of Fa/Fb— 3. Since there is 
no evidence for significant soft X-ray or dust absorption 
(see §2.3, Dai et al. 2006, Poindexter et al. 2008) this dif- 
ference in flux ratios must be due to microlensing. In §3 
we proceed with fitting the light-curves of HE 1104-1805 
with a sophisticated microlensing model. The microlens- 
ing analysis of the light-curves of RX J 1131-1231 will be 
presented in a companion paper by Dai et al. (2009). As 
with RX J1131 — 1231, we see no evidence for variations in 
the spectral hardness, largely due to the large uncertain- 
ties. 

3. X-RAY MICROLENSING IN HE 1104-1805 

We modeled the microlensing in HE 1104—1805 as in 
Morgan et al. (2008a) and Poindexter et al. (2008) us- 
ing a fixed 162 day time delay, the i?-band light curve 



constructed from our SMARTS data (Poindexter et al. 
2007), OGLE (Wyrzykowski et al. 2003), and Ofek & 
Maoz (2003), and the delay-corrected X-ray flux measure- 
ment. The main free parameters of our lensing model for 
HE 1104—1805 are the X-ray and UV (rest-frame) scale 
lengths, a microlcns mass scale, a mass fraction of the 
local surface density comprised of stars, and a velocity 
vector describing the motion of the AGN emission re- 
gion across the microlensing caustics. We used the lens 
model sequence from Poindexter et al. (2007), where we 
start from a constant mass-to-light ratio model for the 
mass distribution, defined by Jm/l = 1, an d then reduce 
its mass in 10% increments while adding an NFW halo, 
where Jm/l = would correspond to a pure halo model. 
The measured time delay requires /m/l — 0.3. We made 
8 random realizations of the star fields near each image 
and then generated 8192 2 pixel magnification maps us- 
ing the methods of Kochanek (2004). The maps had an 
outer scale of 10i?E = 4.3 x 10 16 cm and an inner scale 
of 10i?E/8192 = 4.3 x 10 13 cm. Based on its emission 
line widths, Peng et al. (2006) estimated that the black 
hole mass in HE 1104-1805 is M BH = 2.4 x 1O 9 M , cor- 
responding to a gravitational radius of r g = GMbh/c 2 = 
3.6 x 10 14 cm that is well-resolved in the magnification 
maps. We modeled the surface brightness of the emission 
regions as a face-on, thin disk (Shakura & Sunyaev 1973) 
without the central temperature depression, 
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jum J V10 9 M Q J \r)L E/ 
is the radius at which the disk temperature matches the 
rest wavelength of the observations, kT — h p c/X Test , M 
is the accretion rate, L/Le is the luminosity in units of 
the Eddington luminosity and r\ = L/(Mc 2 ) is the accre- 
tion efficiency. We scale the results assuming a radiative 
efficiency of r\ = 10%. While the available energy of accre- 
tion is larger, the division of that energy between radia- 
tive, advective and kinetic losses is not directly calculable 
at present. Studies of the growth of black holes tend to 
find that the average radiative efficiency of accretion is less 
than 10% (e.g. Shankar et al. 2007), and low radiative ef- 
ficiencies are also needed to reconcile thin disk theory with 
microlensing accretion disk size estimates (e.g., Morgan et 
al. 2007). 

For the UV (rest-frame) emission we can neglect the 
central temperature depression provided i?A roat >> Rin 
where Ri n ~ 2r g is the inner edge of the disk. While this 
thermal accretion disk model is not appropriate for the 
non-thermal X-ray emission, Mortonson et al. (2005) have 
shown that microlensing essentially measures the half-light 
radius of the emission region, so we used the same profile 
for the X-rays to simplify the computations but regard it 
as a measurement of the half-light radius R1/2 — 2AAR\. 
The microlenses were given a power-law mass distribution 
dN/dM oc M -1 - 3 with a dynamic range in mass of a factor 
of 50, roughly matching the Galactic disk mass function 
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of Gould (2000), and we use a uniform prior on the mean 
microlcns mass of 0.1M© < (M) < M©. We modeled 
the probability distribution for the physical velocities as 
in Poindexter et al. (2008) where our projected velocity 
onto the lens plane is 73 km/s, the velocity dispersion of 
the lens is 308 km/s based on the critical radius of the 
mass model, and the lens and source have rms peculiar 
velocities of 153 and 71 km/s respectively. The data were 
then modeled using the Bayesian Monte Carlo method of 
Kochanek (2004). 

The microlensing data constrains the source velocity and 
size in "Einstein" length units of (M) 1 / 2 cm. We can de- 
termine the true sizes using the prior on the microlens 
masses, the prior on the physical velocities or both. Fig- 
ure 7 shows the probability distribution found for (M) 
assuming the prior on the physical velocities. While the 
distribution peaks at somewhat lower mass (~ 0.05M©) 
than the range of our prior (0.1M© to 1.0M©), the dis- 
tribution is so broad (0.005M© to 0.9M© at l-cr) that we 
can view them as statistically consistent. For example, 
none of the results depend significantly on the inclusion 
of the mass prior. Since the mass scales as (M) oc v 2 /v 2 
where v e is the physical velocity from the prior and v e is 
the Einstein-unit velocity estimated from the microlensing 
models, it is difficult to tightly constrain the mean mass 
with a single lens given the unknown peculiar velocities 
of the lens and source even with perfect knowledge of the 
Einstein-unit velocity. 

Figure 8 shows the raw estimates of the source sizes 
f s in Einstein units, where r s = f s (M/M©) 1//2 . Both of 
these are scaled source sizes R\ defined by the emission 
profile presented earlier in this section. The consequence 
of adding the prior on the microlens masses is to drive the 
source sizes f s to be modestly smaller. In essence, a given 
level of microlensing variability can be created by a large 
source moving rapidly (big f s and v e ) or a small source 
moving more slowly (small f s and v e ). But high effective 
velocities correspond to small mass scales given a fixed 
physical velocity, so imposing the mass prior given the 
distribution seen in Figure 7 will favor lower velocities v e 
and smaller sources r s . It will have less effect on the phys- 
ical source size r s = f^M/M©) 1 / 2 because the smaller 
size r s is balanced by the larger mass scale (M/M©) (see 
Kochanek 2004). While Figure 8 truncates the probabil- 
ity distribution for the X-ray source at the pixel scale of 
the magnification patterns, it would not converge to zero 
if we used magnification patterns with smaller pixel scales 
to extend the distribution to smaller source sizes. This 
is a consequence of not having a well-sampled X-ray light 
curve. The source must be small enough to allow the X- 
ray flux ratio to differ from the optical, but still smaller 
sources are disfavored because even more discrepant flux 
ratios become more likely. However, reducing the size still 
further has no effect on the probabilities because it con- 
verges to the likelihood of the data points avoiding precise 
caustic crossings. Unfortunately, this does mean that we 
will obtain robust upper limits on the X-ray source size 
but not lower limits. 

Figure 9 shows the results for the rest-frame UV and 
X-ray source sizes R\ assuming either uniform or loga- 
rithmic priors for the size. We used R\ for both so that 
the sizes are easily compared. If we regard the source 



size at 2000A as measuring a projected area, then R\ 
scales as (cosz)~ 1//2 for inclination angle i of a thin in- 
clined accretion disk. Similarly, the X-ray size in Fig- 
ure 9 should be increased by a factor of 2.44 to con- 
vert from R\ to the half-light radius R\/2- The disk 
size at 2000A (i?-band in the rest frame at z s — 2.32) 
is 15.6 < log(i?200 \J cos(i) / cm) < 16.2 for a logarithmic 
prior on the size and 15.8 < log(i? 2 oo \J cos(i) /cm) < 16.3 
for a linear prior on the size. This agrees well with our 
earlier estimates in Poindexter et al. (2008) and Mor- 
gan et al. (2008a). The X-ray size is less certain and 
more prior-dependent, with an estimated half-light radius 
14.2 < log(i?x.i/2/ cm ) < 15.0 for a logarithmic prior on 
the size and 14.6 < log(i?x.i/2/ cm ) < 15-3 for the linear 
prior on the size. We note that the difference in the sizes 
is entirely driven by the significant differences between the 
X-ray and optical flux ratios at the time of the X-ray ob- 
servation. 

The microlensing data also prefer lens models with low 
stellar mass fractions, as shown in Figure 10. The prob- 
ability distribution peaks at Jm/l — 0.2, in relatively 
good agreement with the value of Jm/l = 0-3 that agrees 
best with the measured time delay and a Hubble constant 
of Hq = 72 km s _1 Mpc -1 . While this preliminary re- 
sult does not lead to a microlensing constraint on Jm/l 
that is tight enough to break the degeneracies between 
the radial mass distribution and the Hubble constant (see 
Kochanek 2002), it is encouraging that there is agreement 
with the time delays given the expected value of H and 
that the microlensing requires a dark matter-dominated 
mass distribution. A detailed microlensing analysis of 
RX J1131 — 1231 that provides similar constraints on the 
structure of the lens and the structure of the quasar is 
presented in Dai et al. (2009). 

Figures 11 and 12 show an example of a trial that fits the 
data well. In this trial the Einstein unit source sizes are 
\ogR\ = 14.2 and 16.2 for the X-ray and UV rest-frame 
sources respectively and the effective velocity corresponds 
to a sensible microlens mass scale. We chose this solution 
mainly because of the curious origin of the initial i?-band 
peak, where what might be interpreted as a single heavily 
smoothed caustic crossing is actually a pair of competing 
caustic crossings in the A and B images, where the A im- 
age dominates because of its shallower magnification gra- 
dients. The X-ray epoch then corresponds to time when 
the source is in a valley between caustics in both images. 
We note that this particular solution probably would not 
reproduce the wavelength dependence of the optical and 
IR flux ratios in HE 1104—1805 discussed and modeled 
Poindexter et al. (2008). 

4. DISCUSSION AND CONCLUSIONS 

While there is a consensus that the origin of the keV 
X-ray continuum emission from quasars is inverse Comp- 
ton scattering (e.g., see the review of Reynolds & Nowak 
2003), the extent and geometric configuration of the emis- 
sion region is less clear. In the standard disk-corona model 
(e.g., Haardt & Maraschi 1991; Merloni 2003), it is pro- 
duced in a hot, extended corona surrounding the disk. In 
the general relativistic MHD simulations of Hirose et al. 
(2004, also Machida & Matsumoto 2003) the region of the 
corona with the highest current densities lies close to the 
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inner edge of the disk. Under the assumption that these 
are also the regions with the highest dissipation rates for 
heating the corona (see Hirose et al. 2004), they will dom- 
inate the X-ray emission and we would expect the X-ray 
source size to track the inner edge of the disk. 

Microlensing of the X-ray source allows us to observa- 
tionally measure the size of the X-ray emission region. We 
start with well-sampled optical light curves that allow us 
to determine the size of the accretion disk. Morgan et al. 
(2008a) found a clear scaling of the disk size with black 

2 /3 

hole mass that is consistent with the M B ' H scaling of sim- 
ple thin disk theory for a fixed Eddington accretion rate 
(Shakura & Sunyaev 1973). We can estimate the disk size 
by two routes other than microlensing. First we can esti- 
mate it from the observed flux at some wavelength. Size 
estimates from the observed flux do not depend on the 
black hole mass but depend only on the assumed temper- 
ature profile and the assumption that the disk radiates 
as a black body at that temperature. Second, thin disk 
theory predicts the scale length through the expression 
we presented in the beginning of §3. Size estimates from 
thin disk theory depend on the black-hole mass, the Ed- 
dington factor and the wavelength. The disk size from 
our microlensing analysis is roughly consistent with thin 
disk theory if the accretion efficiency is relatively low, but 
both the microlensing and theoretical sizes are significantly 
larger than that estimated from the observed optical flux 
of the quasar and the same disk model. This mismatch 
between the sizes estimated from microlensing and from 
the flux was first pointed out by Pooley et al. (2007). One 
solution to this problem may be to use a shallower tem- 
perature profile than the T oc i?~ 3 / 4 of simple thin disk 
theory, and Poindexter et al. (2008) found some evidence 
for a shallower profile in HE 1104—1805 by examining the 
optical/near-IR wavelength dependence of microlensing in 
this system. 

Our X-ray monitoring observations of RX J1131 1231 
and HE 1104—1805 show significant flux variability in all 
images of both quasars. The factor of ~ 20 X-ray flux 
variability in image A of RX J1131— 1231 is interesting not 
only because of the high amplitude but also because it is a 
highly magnified saddle point image. Schechter & Wambs- 
ganss (2002) pointed out that in a dark matter dominated 
galaxy, microlensing is in a regime where images at min- 
ima and saddle points of the time-delay surface behave 
differently. In particular, the saddle point images should 
show significantly higher variability amplitudes, just as we 
observe in RX J1131 — 1231. The relatively large sizes of 
accretion disks in the optical relative to the typical Ein- 
stein radius of the stars seems to have masked much of 
this effect in the optical signatures of microlensing - for 
example, where the X-ray flux ratio Fa/ Fc changed by a 
factor of 4 during our monitoring period, the optical flux 
ratio changed by only a factor of 1.4. The differences be- 
tween the X-ray and optical flux ratios of HE 1104—1805 
also imply the presence of X-ray microlensing. Neither 
system shows evidence for soft X-ray or dust absorption 
that could explain the differences. 

While we measure the source sizes by modeling the data 
in detail, the essence of the measurement is comparing the 
amplitude of the X-ray microlensing variability to that in 
the optical. If the amplitude for the X-ray microlensing 



is larger than in the optical, then the X-ray source must 
be more compact, and this is what we observed in both 
RX J1131-1231 and HE 1104-1805. Given the estimated 
black hole mass of M B h = 2.4 x 10 9 M© (Peng et al. 
2006), our fits to the X-ray and optical light curves us- 
ing the microlensing model described in §2.5 indicate that 
the X-ray emission in HE 1104—1805 comes from a region 
smaller than ~ 6r g and the UV emission region has a size 
of ~ 30r 9 . We note that the source sizes of HE 1104-1805 
were derived using a Bayesian analysis which requires a 
prior. As shown in Figures 8 and 9 our upper limit on the 
X-ray source size is robust to prior selection but not the 
lower limit. Obtaining a more precise X-ray size estimate 
requires a more extensive X-ray light curve for the system. 
This is somewhat challenging for HE 1104—1805 because 
the long time delay of HE 1104-1805 (162 days, Morgan 
et al. 2008a) and the intrinsic variability we observe here 
means that this lens requires observations separated by the 
time delay in order to correct for intrinsic variability. This 
is less of a problem for almost all other candidates for X- 
ray monitoring because they have far shorter time delays. 
For example, in RX J1131 — 1231 we observe some intrin- 
sic variability but it is significant only on time scales much 
longer than the delays between the A/B/C images. Given 
the Morgan et al. (2008b) estimates for PG 1115+080, 
the Dai et al. (2009) results for RX J1131-1231, and our 
present results for HE 1104—1805, the X-ray microlensing 
results suggest that the size of the X-ray emission arises 
from close to the inner edge of the disk. Thus, our esti- 
mated sizes for the X-ray and UV emission regions place 
significant constraints on AGN models that assume ex- 
tended coronal envelopes around accretion disks. Reduc- 
ing the uncertainties in these measurements and measur- 
ing the scaling of the size with the estimated masses of 
the black holes requires longer X-ray time series and more 
systems, but there is clearly no technical barrier to making 
the measurements and extending them to the dependence 
of the size on X-ray energy. 

Combining the data from the time series also allows us 
to search for extended emission from groups or clusters 
near the lens, which can be important for lens models (e.g. 
Williams et al. 2006, Fassnacht et al. 2007). Here we ana- 
lyzed the X-ray data for the foreground cluster found near 
RX J1131-1231 by Morgan et al. (2006) in more detail, 
and set limits on the presence of any other clusters near 
both systems (see §2.4). We also used the combined images 
to search for central or odd images, although this is chal- 
lenging given the resolution of Chandra and the expected 
faintness of the images. While we reached la limits on the 
flux ratios of > 100, significant constraints on the central 
images require limits better by another order of magnitude 
(Keeton 2003). 

We acknowledge financial support from NASA via the 
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TABLE 1 

Log of Observations of Quasars RX J1131-1231 and HE 1104-1805 



Chandra Exposure 

Epoch Observation Observation Time Wa° A<B a Nc a N D a 

Date ID (ks) counts counts counts counts 



Observations of RX J1131-1231 



1 


2004 April 12 




4814 


10. 


425 ± 22 


2950 ± 54 


839 ± 29 


211 ± 15 


2 


2006 March 10 




6913 


4.9 


393 ± 20 


624 ± 25 


204 ± 14 


103 ± 10 


3 


2006 March 15 




6912 


4.4 


381 ± 20 


616 ± 25 


233 ± 15 


93 ± 10 


4 


2006 April 12 




6914 


4.9 


413 ± 20 


507 ± 23 


146 ± 12 


131 ± 12 


5 


2006 November 


10 


6915 


4.8 


3708 ± 61 


1411 ± 38 


367 ± 19 


155 ± 13 


6 


2006 November 


13 


6916 


4.8 


3833 ± 62 


1618 ± 40 


415 ± 20 


115 ± 11 



Observations of HE 1104-1805 



1 


2000 June 10 


375 


47.4 


1024 ± 32 


547 ± 23 


2 


2006 February 16 


6918 


5.0 


38 ± 6 


86 ± 9 


3 


2006 March 15 


6917 


4.6 


40 ± 6 


59 ± 8 


4 


2006 April 09 


6919 


4.9 


33 ± 6 


61 ± 8 


5 


2006 October 31 


6920 


5.0 


60 ± 8 


40 ± 6 


6 


2006 November 08 


6921 


4.9 


52 ± 7 


57 ± 8 



a Background-subtracted source counts for events with energies in the 
0.2-10 keV band. The counts for images B and C for the 2004 April 
12 observation of RX J1131 — 1231 are corrected for pile-up by factors of 
37% and 16%, respectively. The counts for images A and B for the 2006 
November 10 observation of RX J1131 — 1231 are corrected for pile-up by 
factor of 32% and 13%, respectively. The counts for images A and B for 
the 2006 November 13 observation of RX J1131 — 1231 are corrected for 
pile-up by factors of 33% and 17%, respectively. 
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TABLE 2 

OPTICAL AND X-RAY OFFSETS OF RX J1131-1231 AND HE 1104-1805 IMAGES 



Telescope 


c 

Aa("), A6(") 


A B 

Aq("), A«(") Aa("), A5(") 


D 

Aa("), A<5(") 


RX J1131-1231 Image Offsets 


HST 

Chandra 


0,0 
0,0 


0.588±0.003,1.120±0.003 0.618±0.003,2.307±0.003 
0.59±0.02, 1.12±0.02 0.62±0.02, 2.31±0.02 


-2.517±0.003, 1.998±0.003 
-2.52±0.02, 2.00±0.02 


HE 1104-1805 Image Offsets 


HST 

Chandra 




0.0 2.901±0.003,-1.332±0.003 
0,0 2.92±0.02, -1.33±0.02 
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TABLE 3 

RESULTS FROM FITS TO THE Chandra IMAGE SPECTRA OF RX J1131-1231 



Epoch 



Parameter 



Values For 
Image A b 



Values For 
Image B fa 



Values For 
Image C b 



Values For 
Image D b 



2-10 keV Flux (10~ 13 erg s" 1 cm" 
C — statistic/ nbins 



0.2-2 koV Flux (lO -1,3 erg s" 1 cm" 
2-10 koV Flux (10~ 13 org s _1 cm" 
C — statistic/ 'nbins 



0.2-2 kcV Flux (10~ 13 erg cm" 
2-10 kcV Flux (10~ 13 erg s" 1 cm" 
C — statistic/ nbins 

p(x 2 M c 



13 erg s 1 cm 
2-10 kcV Flux (10" 13 erg s" 1 cm" 
C — statistic /nbins 
X 2 /x 

p( x 2 /vy 



1 44+ oos 

-0.08 

1 4 +0 ' 2 

9 7+0-4 
z - 1 -0.4 

556/787 
25/31 
0.76 

1 58+ 015 

J - oo -0.15 

2 n + 0.6 

4 5+ 1 ' 2 

*'°-0.9 

531/838 
33/28 
0.24 

1 61+ 015 

- L - ul -0.15 
o 1+0.8 
°- 1 -0.6 

4 6 + 1 ' 2 

-0.9 

521/838 
26/27 
0.52 



1.53 



+ 0.14 
0.14 
o n + 0.6 
°-"-0.6 

°' i -1.0 

541/838 
27/29 
0.58 



1 4i+0-° 7 

661/787 
105/97 
0.27 



i.47: 



1.63"* 



D + 0.11 

3 -o.n 

7.oti;? 

653/838 
61/52 
0.18 



1.73- 1 



D + 0.11 

3 -0.11 

6.7ti;5 

555/838 
42/48 
0.73 

1 64+ 012 

3 9+ - 7 

5 4+ 12 
J -^-i.o 

548/838 
40/41 
0.52 



7+0.12 
-0.12 
o 1+0.5 
"^-O.B 
o O+0.6 
d - S -0.5 

583/787 
28/41 
0.95 



1.95! 



1.62"* 



,+0.19 
-0.19 
i O + 0.3 
1 -°-0.3 
9 q+0.7 
^- -0.6 

435/838 
14.4/16 
0.57 

, --+0.18 
1 '°'-0.18 
1 s+0.5 
1 -°-0.5 

n +0.9 

425/838 
18.5/17 
0.36 

1 70+ 23 

1 - lu -0.23 

1 9 + 0.4 
1 - A -0.4 
q r+0.6 
i -°-0.5 

333/838 
7.7/10 
0.66 



,+0.21 
-0.20 

9+°- 2 
u.y_ 2 

7+ ' 2 

' -0.2 

368/787 
15.4/15 
0.42 



1.66"* 



- + 0.28 
"-0.27 
n S+ ' 4 

1 1 + - 5 
1 ' 1 -0.4 

300/838 
2.4/6 
0.9 

i CA+0.28 

i.6or 28 

n s+ - 4 

1 9 + 0.6 
A - z -0.5 

292/838 
3.8/6 
0.7 

1 62+ 024 

1 ' DZ -0.24 

1 0+°' 4 

1 c + 0.6 
i - b -0.5 

350/838 
13/9 
0.15 



r 


1 59+ 07 

1 ' J3 -0 .07 


1 09+0. 08 


1 90 +01S 

lau -0 .15 


1 54+ ' 22 


0.2-2 kcV Flux (10~ 13 erg s" 1 cm" 2 ) 


29+ 2 


13+1 


3 4+ - 5 

J ' -0.5 


1 ,+0.3 


2-10 keV Flux (10~ 13 erg s _1 cm" 2 ) 


35+ 3 


nil 


9 S+ - 7 


1 9+°' 6 


C — statistic/ nbins 


771/838 


626/838 


502/838 


402/838 


x 2 /» 


72/59 


66/63 


34/28 


7/11 


p(x 2 /»r 


0.1 


0.37 


0.19 


0.7 



r 


i ci+0.07 
1 - ol -0.07 


1 09+O.O8 
J--»^_0.08 


1 76 +0 - 14 

!■ 'D-0.14 


9 90 + O.29 
■^ S -0.28 


0.2-2 kcV Flux (10~ 13 erg s _1 cm" 2 ) 


30t 2 


15+ 2 


o 9+0.5 
d -^-0.5 


1 4+0-5 


2-10 keV Flux (10~ 13 erg s _1 cm" 2 ) 


36±| 


13± 2 


3 6+ - 9 
3 -°-0.8 


6 +0 ' 3 

U.O_ 2 


C — statistic /nbins 


768/838 


658/838 


489/838 


276/838 


X 2 /" 


65/59 


62/70 


29/33 


9/8 


p(x 2 /vT 


0.29 


0.76 


0.66 


0.36 



a Model 1 consists of a power-law modified by Galactic absorption. 
b All errors are for 90% confidence unless mentioned otherwise with all 
parameters taken to be of interest except absolute normalization. 
c P(x 2 /t / ) is the probability of exceeding \ 2 f° r v degrees of freedom if 
the model is correct. 
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TABLE 4 

RESULTS FROM FITS TO THE Chandra IMAGE SPECTRA OF HE 1104-1805 



Epoch 


Parameter^ 1 






Values For 
Image A b 


Values For 
Image B b 


1 


r 

0.2-2 koV Flux (10~ 14 erg s _1 
2-10 koV Flux (10~ 14 erg s" 1 
C — statistic / nbins 


cm" 
cm" 


- 2 ) 
- 2 ) 


1 59+ 08 

10.6+2 = 
623/790 


i-««±S:" 

4.28±°;S 
3.97t°o-.8 


2 


V 

0.2-2 koV Flux (10~ 14 erg s" 1 
2-10 keV Flux (10~ 14 erg s" 1 
C — statistic/ nbins 


cm 
cm" 


- 2 ) 
- 2 ) 


1 4 + - 4 

2.6±i:I 

192/838 


10 2+ 6 - 
284/1060 


3 


r 

0.2-2 keV Flux (10~ 14 erg s~ 4 
2-10 keV Flux (10~ 14 erg s _1 
C — statistic/ nbins 


cm" 
cm" 


- 2 ) 
2 ) 


, „+0.40 
1 -°°-0.41 

3 4+ 2 ' 2 

J ^-2.3 

5 9+ 4 '° 

O-y -3.0 

186/838 


2 4+ - 4 
SAtli 

195/1060 


4 


r 

0.2-2 keV Flux (10~ 14 erg s~ 4 
2-10 keV Flux (10~ 14 erg s _1 
C — ' statz sttc / ' nbzn s 


cm 
cm" 


- 2 ) 

- 2 ) 


1 76+ 050 
2 9+ 2 ' 3 

o 9 + 4.0 
°' — 2.0 
1 60 /838 


6-5+1? 

4 Q + 3-5 

224/1060 


5 


r 

0.2-2 keV Flux (10" 14 erg s" 1 
2-10 keV Flux (10" 14 erg s" 1 
C — statistic /nbins 


cm" 
cm" 


- 2 ) 
- 2 ) 


1 50 + 035 
l.OM_ 34 

4.5i 2 - 6 3 
6.9+3° 
236/838 


„ ,0 + 0.49 

z.oo_ 46 

6-0131 
1.8+ 2 ° 
165/1060 


6 


r 

0.2-2 keV Flux (10" 14 erg s" 1 
2-10 keV Flux (10" 14 erg s" 1 
C — statistic /nbins 


cm 
cm" 


- 2 ) 
- 2 ) 


1 61+ - 36 
44 +2.6 

6 4+ 4 ' 7 
230/838 


i-8i±°o: 4 3? 
4.6i 2 ; 3 

4 7+3-6 
' -2.4 

222/1060 



a Model 1 consists of a power-law modified by Galactic absorption. The 
X-ray fluxes represent unabsorbed values. 

6 All errors are for 90% confidence unless mentioned otherwise with all 
parameters taken to be of interest except absolute normalization. 
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Fig. 1. — Lucy-Richardson deconvolved images in the 0.2 - 10 keV bandpass for the Chandra observations of RX J1131— 1231 (left) and 
HE 1104—1805 (right). The brightness scale of the images is proportional to the count rate. The images are displayed with a linear brightness 
scale. East is to the left and North is up. 
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Fig. 2. — Adaptively smoothed image of the 10 ks 2004 Chandra observation of RX J1131— 1231. To reduce background contamination 
and to enhance possible soft extended X-ray emission we filtered the image to include only photons with energies ranging between 0.4 and 
3.0 keV. 
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Fig. 3. — Stacked image of the six observations of RX J1131— 1231 listed in Table 1. The stacked image was filtered to include only photons 
with energies ranging between 0.4 and 3.0 keV and adaptively smoothed. We have indicated with a cross the location of the reported, but low 
significance, extended X-ray emission claimed to possibly originate from a lensing cluster of galaxies. There is no indication of any extended 



emission. 
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Fig. 4. — Stacked image of the six observations of HE 1104—1805 listed in Table 1. The stacked image was filtered to include only photons 
with energies ranging between 0.4 and 3.0 keV and adaptively smoothed. There is no indication of any extended emission originating from a 
group or cluster of galaxies. 
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0.001 1 — 1 — 1 — 1 — 1 — ' — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 
3000 3500 4000 4500 

HJD - 2450000 (days) 

Fig. 5. — The 0.2—10 keV (large point with error bars) and _R-band (small points) light curves corrected for the time delay of (a) images A, 
B, C and D of RX J1131-1231 and (b) images A and B of HE 1104-1805, respectively. The optical fluxes of RX J1131-1231 are normalized 
such that the normalized optical flux of image D is equal to the 0.2-10 keV count-rate of image D near day 3000. The optical fluxes of 
HE 1104—1805 are normalized such that the normalized optical flux of image A is equal to the 0.2-10 keV count-rate of image A near day 
3780. 
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2000 3000 4000 

HJD - 2450000 (days) 

Fig. 6.— (a) X-ray (0.2-10 keV, large points) and optical (i?-band, small points) flux ratios F A /F C , F B /F C and F D /F C for RX J1131-1231. 
(b) X-ray (0.2-10 keV) and optical (fl-band) flux ratios F A /F B for HE 1104-1805. The optical flux ratios of RX J1131-1231 and 
HE 1104-1805 are corrected for the time-delays. The 0.2-10 keV time-delay corrected flux ratios of RX J1131-1231 and HE 1104-1805 are 
shown with solid triangles and the observed ratios are shown with open triangles. 
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0.0001 0.001 0.01 0.1 1 10 



(M) (M ) 

Fig. 7. — Probability density of the microlens mass scale (M) given our model for the estimated physical velocities. The peak is consistent 
with our alternate prior on the mass scale (O.IMq < (M) < Mq) but broader because it is difficult to tightly constrain the mass scale using 
the velocities. 
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10 14 10 15 10 16 10 17 

r s (<M/M > 1 / 2 cm) 

Fig. 8. — The probability density of the scaled source sizes R\ in Einstein units of (M/Mq) 1 / 2 cm. Here the pixel scale and the mean 
Einstein radius correspond directly to the magnification patterns. The solid curves use only the velocity prior while the dashed lines add the 
mass prior. The mass prior eliminates low mass, high velocity, large source size solutions. Given our single X-ray data point, the probability 
distribution for the X-ray source would converge to a non-zero constant if we extended the calculations to still smaller sources using higher 
resolution magnification patterns. 
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R x (cm) 

Fig. 9. — Disk size at 2000A (top) and X-ray (bottom) for HE 1104—1805, where the solid (dashed) curves assume logarithmic (linear) 
priors on the source size. We show the values of R\ for the assumed surface brightness profile so that the size ratios will be correct. Physically, 
the disk size at 2000A can be increased by a factor of (cosi) 1//2 for the inclination of the disk i, and the X-ray size should be regarded as 
a measurement of the half-light radius with R1/2 = 2.44 lt\. Vertical lines mark the gravitational radius r g = GMbh /c 2 and accretion disk 
size given the wavelength, the estimated black hole mass, and Eddington limited accretion with an efficiency of r\ = 0.1. The final vertical 
line marks the Einstein radius for solar mass stars. 
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fraction of constant M/L f M/L 

Fig. 10. — Constraints on the mass distribution. The light solid, dashed and heavy solid curves show the constraints on the mass of the 
stellar component of the lens mass model based on the microlensing models (heavy solid), the time delay measurements (dashed) and the 
combination of the two (light solid). A constant M/L model has /m/l = 1 antl a pure dark matter halo model has Jm/l = 0j with the stellar 
mass reduced from the constant M/L model in proportion to /m/l 
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1 2 3456789 

Fig. 11. — Simulated microlensing caustics for HE 1104—1805. These show small sections of one pair of the Jm/l = 0.3 magnification 
patterns for images A (left) and B (right) where higher magnifications are white and lower magnifications are black. The lines show the 
path followed by the source for the period covered by the ij-band data starting from the circled point. The large circle is the size of the UV 
(rest-frame) source (log R\ = 16.2) and the small (almost invisible) circle is the size of the X-ray source (logR\ = 14.2). 
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Fig. 12. — The ij-band and X-ray magnitude differences for the trial shown in Fig. 11 as compared to the observations. This trial was 
chosen because it fit the data well and we found the origin of the initial peak as a pair of competing caustic crossings in the two images 
intriguing. The dramatic differences in the vertical scales show that an X-ray light curve of this lens should be much more dramatic than the 
ij-band. 



